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We discuss details of the preparation and detection of entangled electron-nuclear spin states in 
15 N@C6o together with a quantitative evaluation of the complete density matrix. All four Bell 
states of a two qubit subsystem were analyzed. In addition we find a quantum critical temperature 
of T q = 7.73 K for this system at an electron spin resonance frequency of 95 GHz. 
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I. INTRODUCTION 

In recent years quantum information theory has come 
to exciting new ideas in the fields of quantum teleporta- 
tion 0, Q .quantum cryptography Q and quantum com- 
putation |iJ. l| . It was stimulated enormously by the 
discovery of the factoring |(| Q and searching algorithms 
0, which demonstrate that a quantum computer is 
capable of tasks, that are impossible to solve on a clas- 
sical computer. Consequently, over the past few years 
various physical systems were examined or proposed for 
their potential use as hardware for quantum computing. 
These include ions confined to an electromagnetic trap 
[HI EH photons 0] > doped semiconductors 0], Cooper 
pair states in superconductors 0,0, Josephson junc- 
tions 16, , and nuclear spins in liquids (NMR quantum 
computing) fll ITl IM EHl I2I I2H li^fl . 

One of the most interesting components in many quan- 
tum algorithms are nonlocal quantum correlations that 
violate our conception of the classical world. These so 
called entan gled states have been discussed for a long 
time [M Hi |H 0. Experimental preparations of 
entangled states of photons have been published |30l l31j . 
The pseudo entanglement of three coupled nuclear spins 
in the so-called Greenberger-Horne-Zeilinger (GHZ) 0] 
state was demonstrated by NMR [33. 1.3 ll 13a . l3a |. En- 
tanglement was also achieved among ions confined in an 
electromagnetic trap 0. 

In this contribution we present procedures for creat- 
ing and detecting entanglement between electronic and 
nuclear spin states in the solid phase of the endohcdral 
fullerene 14 N@C6o which has already been proposed as 
a basic unit for a scalable quantum computer |38l l39| . 
We will show, however, that 15 N@Ceo is not a qubit 
but rather a multi-qubit system and specific addressing 
schemes must be applied when considering all quantum 
levels. From its multiple quantum states we will project 
to a two qubit subsystem and demonstrate the generation 
of entanglement within this subsystem. The experimen- 
tal app roach is related to the preliminary work presented 
in [4(j where entanglement has been simulated between 



an electron spin 1/2 and a nuclear spin 1/2 of a radical 
electron in a single crystal of malonic. Due to the strong 
hyperfine interaction in the system presented here entan- 
glement can be achieved by a factor of thousand shorter 
time scale than in liquid state NMR. A brief account of 
our approach to 14 N@C6o was already published in [4lj . 
Here we present further details on the preparation and 
tomography. 

The experiments reported here are carried out in a 
mixed spin state of an ensemble of spins. Since the com- 
plete density matrix in such a system is separable strictly 
speaking entanglement like in pure quantum states can 
not be achieved |42l liflj . However, the unitary opera- 
tions applied to the mixed states density matrices are 
identical to those applied to pure systems and moreover, 
the resulting density matrices have the same operator 
structure as the ones of the pure system except for a dif- 
ferent overall factor. Therefore these states are called 
pseudo pure or pseudo entangled [l^, Hfl |21j . In the end 
we will show that these restrictions can be overcome for 
the system under investigation for high enough magnetic 
fields (95 GHz electron spin resonance frequency)and low 
enough temperature( T = 7.73 K) with current technol- 
ogy. The quantum limit where pseudo entangled states 
become entangled states is therefore well in reach. 

Before we start to present the procedures for the more 
complicated case of 15 N@C6o we sketch as a reminder 
how entanglement could be achieved in a system of two 
coupled quantum bits (qubits) represented for example 
by an electron spin S — 1/2 and a nuclear spin 1=1/2. 
Starting from the Zeeman product states 



\m s mi) = I TT), I tl), I IT), I U) , 



(1) 



entanglement is achieved by applying first a Hadamard 
transformation H on one spin followed by a controlled 
not (CNOT) operation. Suppose we start from the pure 
state I |t) we can create the entangled state \E' + by 
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H J_ 

~* y/2 



II) + 1 IT) 



CNOT J_ 



I Tl) + I IT) • (2) 



Depending on the applied unitary transformations and 
initial states, all four entangled states (Bell-basis) of a 
two qubit system can be generated by this procedure: 



'Electronic address: m.mehring@physik.uni-stuttgart.de 



1 



1 



= - (|TT> ±IU», |* ± ) = -(|TI)±|IT)) ■ (3) 
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As we will see later in detail those transformations can 
be realized by selective high frequency pulses applied to 
allowed transitions of the system. We will use in the fol- 
lowing the qubit notation \ms mi) = 1 00) , 1 01) , 1 10) , 1 1 1) 
as well as the state labelling |1), |2), |3), |4). 



a free nitrogen j47j | atom because of the confinement in- 
side the carbon cage. For strong enough magnetic fields, 
where \a\, luj -C wj, second order contributions to the hy- 
pcrfine coupling can be neglected and the Hamiltonian is 
in first order given by: 
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FIG. 1: The eight quantum states of a two spin system with 
S — 3/2 and I — 1/2. The arrows denote allowed transitions. 
The dotted arrows indicate forbidden transitions that belong 
to the entangled states discussed in this publication. The 
separation of the energy levels is not to scale. 

The experiments were performed on a powder sample 
of 15 N@C 6 o- The 15 N atom resides in the center of the 
Ceo molecule. There is no charge transfer from the atom 
to the cage. The nitrogen atom is paramagnetic due to 
the half filled p orbital shell with three unpaired electrons 
which form a total electron spin S = 3/2. Because of the 
15 N isotope enrichment the nuclear spin of the nitrogen 
atom is / = 1/2. The magnetic resonance properties of 
the sample have already been discussed in detail in 
Irsj . The nitrogen atoms were implanted into the carbon 
cage by Weidinger and co-workers by simultaneous 
evaporation of Ceo and ion bombardment onto a cooled 
target. The sample was purified by high-pressure liquid 
chromatography. 

With the magnetic field oriented along the z-axis, the 
Hamiltonian of this system becomes 



H = h (u>sS z + loiI z + 27raSI) 



(4) 



where u>s and oji are the Larmor frequencies of electron 
spin S = 3/2 and the nuclear spin / = 1/2. Due to the 
high symmetry of the molecule the hyperfine interaction 
a = —22.08 MHz is isotropic. This value is larger than for 



H = h (ujsS z + ujiI z + 2na S z I z ) 



The eight eigenvalues of this system are 



(5) 



En 



h(ivsrns + uimi + 2irams mi) , (6) 



depending on the spin quantum numbers m§ — ±1/2, 
±3/2 of the electron spin and mi = ±1/2 of the nuclear 
spin. The eight eigenstates of the system are given by 



\m 8 mj) = |1),|2),...,|8) 
3 1 



(7) 



The corresponding energy level scheme is shown in Fig. 
n Allowed transitions are indicated by arrows. There are 
in first order three degenerate electron spin transitions 
(Ams = ±1) for mi = +1/2 and and another three for 
mi = —1/2. Because of this degeneracy the transitions 
can only be excited simultaneously. There are, however, 
four nuclear spin transitions (Ami = ±1)- The resulting 
electron spin resonance (ESR) and electron nuclear dou- 
ble resonance (ENDOR) spectra shown in Fig. [21 verify 
this level scheme. The lines are labelled according to the 
spin states mi and ms- The ESR spectrum exhibits two 
lines separated by 22.08 MHz. The ENDOR spectrum 
consists of four lines located at 9.67 MHz, 12.41 MHz, 
31.69 MHz and 34.54 MHz. 



B. Entangled spin states in 15 N@C 6 o 

Instead of the four entangled Bell states for two cou- 
pled spins 1/2 the endohedral fullerene 15 N@C6o allows 
for 2 x 12 entangled states of the type 
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71 



(li)±|fc» 



(8) 



with jk £ {14,23,36,45,58,67,16,25,38,47,18,27}. 
Here we will restrict ourselves to the preparation and 
detection of the following two examples 
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FIG. 2: X-band ESR (top) and pulsed-ENDOR spectra 
[4^ (bottom) of a powder sample of 15 N@Ceo in a magnetic 
field of 340 mT. The ESR spectrum was obtained by Fourier 
transformation of the electron free induction decay (FID). 
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(9) 



± 



(10) 



= (|00) ± |11)) 

Except for the electron spin state with ms = ±3/2 
these states are equivalent to the two qubit Bell states of 
Eq. ©• The corresponding coherent superpositions are 
indicated in Fig. ^ by dashed arrows. This defines the 
fictitious two state subsystem of the electron spin with 
| ± 3/2) as one qubit. For the second qubit we consider 
the | ± 1/2) states of the nuclear spin. In other words the 
quantum states | ± 3/2 ± 1/2) represent the four states 
of our two qubit subsystem. Of course there are other 
combinations of quantum states that could be regarded 
as an independent two qubit subsystem. In this sense 
the quantum system of 15 N@Ceo represents a multi qubit 



system. In the following we will use the qubit notation 
1 00), |01), 1 10) and |11) for the states | + |+|) lN 



| + \) and | - | - i) 



C. Z-rotations of the quantum states 

In the tomography of quantum states we will apply ex- 
tensively phase rotations about the z-axis. Therefore we 
investigate here briefly how the different quantum states 
behave under these phase rotations. A rotation of a spin 
state of spin 5 by angle <pi and spin / by phase </?2 around 
the quantization axis (z-axis), corresponds to the unitary 
transformations 

U Sz = e~ lViS * and U Iz = e^ 2 ** . (11) 

Under these transformations the quantum state \ms mi) 
exhibits the following phase variation: 



U s , Uijmsmi) = e-^'e-^^'lmsmj) 

= e -^ msVl+mi ^\m s mi) 



(12) 



Let us apply these transformations to a superposition 
state of qubit 1 (subsystem with mj = ±3/2) which re- 
sults in 

(|00) + |10>) (<00| + (10|) 100)(00| + |10)(10| 

+ e- j3vi |00)(10| + e i3v>1 |10)(00|. (13) 

Note that in the case of an electron spin 5=1/2 the 
factor 3 in front of (pi should be replaced by 1 40] . Con- 
sequently the superposition state of nuclear spin states 
with m/ = ±1/2 transforms like 



(|10> + |11»«10| + <11|) 
+ e 



-rot. 



|10)(10| 



|H)(H| 



;-^ 2 |10)(ll| +e 4V2 |ll)(10|. (14) 

Note, however, that applying phase rotations to entan- 
gled states leads to a different behavior as is demon- 
strated here for the Vf^ aim ^18 states: 

(|oi) + |io» ((oi| + (io|) z --^- 101><011 + |10)(10| 
+ e -<(3*'i-*'0|oi)(io|+e^ 3 *' 1 - (pa )|10)(01| (15a) 



(|00) + |11)) ((00| + (11|) 

+ e -*(3<Pi+*>a)|00)(ll| 



-rot 



|00)(00| + |11)(11| 
e <(3*>i+!pa)|ii)(oo|. (15b) 



The entanglement is evidenced by the combined phase 
dependence where both phase angles appear as a sum or 
difference depending on the type of entangled state. In 
contrast the superposition states of individual spins show 
only their corresponding single phase behavior. This will 
be utilized to distinguish different quantum states in the 
following. 



4 



II. EXPERIMENTAL 

The experiments were performed on a diluted pow- 
der sample of 15 N@C6o with a home-built X-band pulsed 
spectrometer including a slotted tube resonator at about 
9.5 GHz. A radio frequency (rf) coil for the excitation of 
the ENDOR transitions was inserted in the homogeneous 
region of the microwave held with frequencies ranging 
from — 40 MHz. Both, the microwave channel at X- 
band frequencies as well as the rf-channel were equipped 
with arbitrary waveform generators (AWG) which al- 
lowed for arbitrary pulse and phase modulation. The 
quantum states were controlled by applying special uni- 
tary transformations implemented as transition selective 
microwave and radio frequency pulses with rotation an- 
gle (3. These pulses are labelled fj-y (/?) for the two elec- 
tron spin transitions corresponding to mj = ±1/2 and 

Px-y (/3) for radio frequency pulses on nuclear spin tran- 
sition between levels j and k (see Fig. QJ. The subscripts 
x or y indicate the phases of the pulses in the correspond- 
ing rotating frame. The selective pulses can be written 
in terms of fictitious spin 1/2 operators FLJ or F^fy', 
which belong to the transition (±) or (jk) |49l l5(| as 



Pulses with arbitrary phase are expressed as 



_ - t (3F± 



(16) 



(17) 



P^ 2 (P) = e-^^'e-^^'eW" . (18) 

All experiments in this contribution were performed at 
a temperature of T = 50 K. 



-ifoFi" p -i0Fi k jfaFi* 



py/kt,! ; At, 


y 

_ B 


i 


ft 


detection 
A« s |f|k™i 














t 








y 

r-n 






I t 
i 


preparation 


entangl. 




! t ' 

tomography i detection 



FIG. 3: Pulse sequence for preparation and detection of the 
pseudo entangled state in 15 N@C60- The left part refers 
to the preparation of the pseudo pure state. The middle 
part represents the preparation and detection of the entan- 
gled state. The FID after the selective ESR pulse (right hand 
side) serves as a monitor signal. The angles ifii and 1^2 of the 
tomography sequence are varied in increments (see text) and 
are used for separating signals of entangled states from those 
of individual spins. The waiting times between the pulses 
were chosen Ati — 5 ps, At2 = 100 fjs, At3 — 150 /is and 
t = 40 ns (see text). 

These pulses serve different purposes as is sketched in 
the general pulse sequence displayed in Fig. [31 The to- 
tal sequence is separated into four segments, namely (1) 



preparation of pseudo pure states, (2) creation of entan- 
glement, (3) state tomography and (4) detection via the 
electron spin free induction decay (FID). 



III. PSEUDO PURE STATES. 

A. Preparation of pseudo pure states 

Like in liquid state NMR quantum computing [l8L I20L 
l2l| we are dealing here with mixed quantum states and 
moreover we start initially from a thermal equilibrium 
state, namely the Boltzmann state which can be ex- 
pressed as 

? B = Tr{e-/W-} (19) 

with /3b = h/k^T and where we have applied the 
approximation 3> wj,<z. With all experiments per- 
formed at T = 50 K and hu s /(k B T) w 0.01 < 1 the 
expression can be further simplified by applying the high 
temperature approximation which results for our S = 3/2 
and 1=1/2 system in 

PB w i (I 8 - K B S Z ) with K B = ^. (20) 

8 KB 1 

and where Is represents the 8x8 identity matrix. Eq. Ij20|l 
can be rewritten as 



PB 



' '1 - \k b ! is 



1 



-KbPb with 



Pp - 7 7 8 _ 2 Sz 



(21) 
(22) 



This defines the 8x8 so-called pseudo Boltzmann ma- 
trix pp which will be the starting point for preparing 
pseudo pure initial density matrices. 

We treat in the following the preparation of the pseudo 
pure density matrix ppio, from which we will extract 
finally the two qubit density matrix pio with diagonal 
components {0,0,1,0}, as an example. It is important 
to prepare a pseudo pure state which not only resembles 
the corresponding pure quantum state with the same op- 
erator structure but also corresponds to a large nuclear 
spin alignment which is a prerequisite for reaching a high 
degree of entanglement. According to Fig. |21we first ap- 
ply a pulse Pt\Po) with fa = arccos(-l/3) = 109,47° 
at the electron spin transitions for mi = —1/2. After a 
waiting time of n = 5 (is all transverse components have 
decayed and a pulse P y ' (ao) pulse with a = tt/2 at the 
1 <-> 2 ENDOR transition follows in order to equalize the 
populations of level 1 and 2. After an additional waiting 
time of T2 = 100 ps all transient components of the den- 
sity matrix have decayed and the pseudo pure state ppio 
has been reached with 



in-ui = {0,0, o,l, i,|,l,0 



(23) 
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where here and in the following we express diagonal 
density matrices by the diagonal element vector where 
the elements in bold type correspond to the matrix pio 
of the subsystem mj = ±3/2 and mi = ±1/2. 

The pscudopure state p\\ can be prepared in a similar 
way by exciting the electron spin transition with m/ = 
+1/2. 




0.5 1.0 1.5 

time (us) 

FIG. 4: Rabi oscillations at the electron spin transitions with 
mi = —1/2. Top: Boltzmann state pp. Bottom: After appli- 
cation of a Py \Po) pulse. 



B. Tomography of pseudo pure states 

The reconstruction of the density matrix of a quan- 
tum state by density matrix tomography has already 
been applied in liquid state NMR quantum computing 
EJ I2H l33t l36l 0]). Here the situation is different be- 
cause we are dealing with ESR and ENDOR transitions 
in a solid at low temperatures. In order to verify the 
proper preparation of the pseudo pure density matrices 
we have applied Rabi oscillation measurements on dif- 
ferent electron and nuclear transitions. This is shown 
for the pseudopure state pio in Fig. 0] as an example. 
There the magnitude of the electron spin free induction 
decay (FID) signal is plotted after applying pulses of vari- 
able length at particular ESR or ENDOR transitions. 
The Rabi oscillations decay due to the microwave- or rf- 
field inhomogeneity. However, their amplitude and initial 
phase represent the population difference in magnitude 
and sign of the particular transition where the pulse is 
applied. The amplitude of the Rabi oscillations was de- 
termined by Fourier transformation and integration over 
the line in the corresponding spectrum. The experimen- 
tal values were calibrated with respect to the Boltzmann 
state pp whose amplitude is defined as Ap = 1/2. Similar 
measurements were obtained for state p\\. 



Figure 0] shows Rabi oscillations at the ESR transition 
with mi = —1/2. A comparison is made between the 
Boltzmann state (top) and an inverted state (bottom) 
after the application of a Py~^ (/3o) pulse. Note the change 
in amplitude and sign with respect to the Boltzmann 
state. The level of inversion was determined to -0.329 
which results in an ang le fa = arccos(-0.329) = 109.2° 
to be compared with the ideal value (3q = arccos(— 1/3) = 
109.5°. 

In addition we performed similar Rabi precession ex- 
periments at the ENDOR transitions 1 *-* 2 and 7^8 
of the pseudo pure state ppio- From the amplitude of the 
Rabi oscillations we have determined the corresponding 
population differences. The deviation from the theoreti- 
cally expected values can be expressed in terms of devia- 
tions of the angles (3o and ao from their ideal values. By 
comparing the expected signal strength 



A pio = -^cosa (l 



cos/? ) 



(24) 



(12) 

with and without the Py 



(ao) pulse leads to an ef- 
fective value of ao = 1.55 = 88.8° which reflects a rather 
small deviation from the ideal value ao = tt/2 = 90.0°. 
Deviations of the angles ao and (3o from their ideal val- 
ues are not due to misadjusting but are rather caused 
by microwave- and rf-field inhomogeneities. This affects 
the diagonal elements of Pp^ and leads to the following 
values 



(exp) 

Pvio 



-0.01,0.01,0, 0.33, ^, 0.17, 1,0.00 \ (25) 



where negative values arise from experimental errors. 
After such a data analysis we obtain for the pseudo pure 
fictitious two qubit density matrix pio 



Pw 



-0.01 

0.01 

10 

0.00 



(26) 



The errors are in the last digit. Note that the num- 
bers 1 and \ derive from pp of the Boltzmann state (see 
Ea. (|21|) ') and are not changed by the preparation pulses. 
Matrix elements of those types are written in the follow- 
ing as bare numbers or fractions. In a similar way as the 
pseudo pure density matrix pn was obtained. 

These pseudo pure states serve as initial states for the 
generation of entangled states, pio was used to generate 
the Bell-states and pu to generate $f s . 



IV. PSEUDO ENTANGLED STATES 

A. Preparation of pseudo entangled states 

Pseudo entangled states were prepared according to 
the pulse sequence depicted in Fig. [3] (central part). In- 
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stead of the usual Hadamard transformation a selective 
it/2 pulse on a specific nuclear spin transition was ap- 
plied which has a related effect as the Hadamard trans- 
form. The CNOT-gate was implemented by a selective 
7r-pulse on an electron spin transition. Either vp* or <I> ± 
states were generated depending on the pseudo pure ini- 
tial state and the transitions used. In order to distinguish 
the different Bell states we have used phase rotations of 
the pulses in the tomography sequence (see Fig. [3J) as 
presented before |io| . 

According to equations © and (fTUf) we need to prepare 
the following density matrices 



with 



f 



p% 7 ± = -\2±7) (2±7| 



and 



P<t>± 



|1±8) (1±8| . 

(27) 

Starting from the pseudo pure state ppio first a ot\ = 
7r/2-pulse was applied at the nuclear spin transition 7 <-> 
8 followed immediately by a f3\ — 7r-pulse at electron spin 
transitions with mj = — 1/2 (see Fig. |3J) leading to the 
following sequence of unitary transformations 

Pp*± = C±Vio (Of) f with Ul 7 = P- (ttt) Pj 8 (+7T/2 

(28) 

Under ideal conditions this would result in 
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(29) 



where the elements in bold type again belong to the 
fictitious two qubit submatrix 
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(30) 
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(32) 



where the elements in bold type again belong to the fic- 
titious two qubit submatrix 
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o ±\\ 




o o h 



(33) 



Here we have used a pulse length of 1.6 fj,s for the ir/2 
radio frequency pulse and 88 ns for the microwave tt- 
pulse. By this entanglement was prepared within ap- 
proximately 1.7/zs starting from the pseudo pure state. 

Note that the density matrices p^± (Eq. $60\ ) and 
•p$± (Eq. correspond to the Bell states defined in 

Eq. {("J) . Here we investigated which density matrices are 
obtained under ideal conditions. In the next section we 
analyze the experimentally obtained density matrices by 
performing a density matrix tomography. 



B. Density matrix tomography 

In the density matrix tomography used here we com- 
bine Rabi precession to determine the populations of the 
quantum levels together with phase incrementation mea- 
surements to discriminate between different off-diagonal 
components. 

Since entangled states are not directly observable we 
need to transform them to an observable state. We there- 
fore apply an entangled state detector consisting for the 
states \l/ ± the pulse sequence Px (— (3%) with 0% = ir fol- 
lowed by Px 78 \a) with a = ir/2 (see Fig. |2J. We note 
that we always use the same angle (3± in the preparation 
and the detection sequence or technically speaking we 
use the same amplitude and pulscwidth. 



The two qubit submatrix of the density matrix p 2 J <s ,± 
is almost identical to p^ ± (see eq. $27 f ) except for addi- 
tional diagonal elements outside the fictitious two qubit 
subsystem. 

Similarly we have prepared the p]^± states by starting 
from the pseudo pure state ppn and applied a 7r/2-pulse 
on nuclear transition 7 <-> 8 followed immediately by a 
7r-Pulse on electron spin transitions with m; = +1/2 

p P 8 $± = t4Vn (U 1 ^ with Ui 8 = P+(±tt) PJ 8 (+tt/2) 

(31) 



1. Phase rotation of entangled states 

An essential aspect of our type of tomography for dis- 
tinguishing between the different Bell states is their de- 
pendence on z-rotations as discussed in section 11*01 40] . 
Accordingly we have implemented the detection sequence 
with P 2 r(— 7r, <j)\) and P 78 (ir/2, fa) pulses. The cor- 
responding unitary transformation of the detection se- 
quence then reads 

Uffa, 2 )=PJ 8 (tt/2, <h)P*(-«, <f>i). (34) 



The measured quantity after the detection is the z- 
magnetization of the electron spin subsystem with mi — 
— 1/2 in the case of ^f T . The observed signal strength is 
therefore expected to vary as 



2) = - 



Tr{F-Ufp^ ± Ul 
Tr{(F z ') 2 } 



27t\ 



l±lcos(3^-0 2 ) 



(35) 
(36) 



where F~ is the fictitious electron spin 3/2 of the m/ = 
— 1/2 subsystem. The phase dependence with (30i — 02 ) 
is characteristic for the \t ± states. The variation with 
30i originates from the ±3/2 levels. 
Similarly we have used the sequence 



Ui 



(37) 



for the detection of the states resulting in a detector 
signal 



, 1 



Tr { F +U™p% ± uf} 



Tr 



{(^ + ) 2 } 



2 3 

tj-g ± 20 cos ( 3 0i + <fo) 



(38) 
(39) 



for and where preparation and detection pulses were 
fT/2 ESR transitions. The phase 



is characteristic for the <f> 



IN 



applied at the mj = 
variation with (30i 
states. 

The phase dependence of the entangled state was mea- 
sured by repeating the whole pulse sequence shown in 
Fig. |3 for different values of cpx and (f2- The angles 
ip 1 and if2 were incremented in steps of A(pi and Atp2, 
leading to phase angles 



tfi(n) — n Aipi — nv\ At 
P2 (n) — 71 Aif2 — nv2 At 



(40) 
(41) 



after n steps. We have introduced here a virtual time 
scale t = nAt which defines the virtual frequencies V\ 
and V2- By incrementing both phase angles 0x and 02 
simultaneously an oscillatory behavior of the detection 
signal according to eqns. I|36|) and l|39|l is expected. In 
the following we call the corresponding traces, shown in 
Fig-© phase interferograms. 

Figure © shows phase interferograms and the corre- 
sponding Fourier spectra for phase incrementation the 
two frequencies 1^1 = 2.5 MHz and V2 = 1.0 MHz. The 
presence of entangled states is revealed by dominant lines 
at 3^1—^2 for the Pp^± states ((c) and (d)) and at 3^i+i/2 
for the Pp^± states ((e) and (f)) in contrast to the single 
spin phase variations seen in (a) and (b). 

The additional weak lines appearing at v\ and 2v\ =p v\ 
are artifacts due to pulse imperfections. A theoretical 
analysis shows that the ratio of the amplitudes of the 




0.5 1.0 1.5 
time (jxs' 



5 10 
frequency (MHz) 



FIG. 5: Left: Phase interferograms for a complete set of Bell 
states according to the pulse sequence in FigEJwith phase in- 
crementation frequencies v\ = 2.5 MHz and v\ = 1.0 MHz. 
The experimental traces belong to pp'Jp- (a) to (c), pp^, + (d), 
pp^+ (e) and pp* 5 - (f). In the first two interferograms only 
one of phase angle was incremented: (a): incrementation of 
02, (b): incrementation of 0i. In the interferograms (c) to 
(f) both phase angles 0i an d 02 were incremented simulta- 
neously. Right: Corresponding power spectra after Fourier 
transformation . There is a dominant line at 3vi — U2 for the 
P P 7 *± states and at 3^i + V2 for the pp 8 $ ± states. See text for 
additional weak lines. 



lines at 2u\ =p v\ and Zv\ =p v\ depends on the deviation 
e/3 of the electron spin rotation angle (3\ from 7r 



A{2v 1 T- V2) 
A(3i/i T v 2 ) 



24 



(l6 - 21 e 2 



with 



Pi 



(42) 

where terms up to second order in ep were considered. 
This we have used to quantify the precision of the applied 
electron spin pulse (3\. For example the mean deviation 
from the ideal situation for state /9p^_ (Fig. 0(c)) can be 
determined to be je^l = 0.23. This is in accordance with 
the fact that a perfect 7r-pulse does not exist, because 
microwave field inhomogeneities alone already prevent a 
complete inversion. ep therefore comprises positive as 
well as negative deviations. Smaller values of ep could 
therefore be achieved by more homogeneous microwave 
fields. The determination of ep as a measure of the pre- 
cision of the entangled states was performed also for the 
other Bell states and is an essential part for the estima- 
tion of the experimental precision of the created states. 
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2. The diagonal part 

The diagonal part of the density matrix of the entan- 
gled states was determined by Rabi oscillations in a sim- 
ilar way as already discussed for the pseudo pure states 
in section IIIII Here we combine electron spin Rabi os- 
cillations together with nuclear spin Rabi oscillations at 
the relevant transitions 1 <-> 2 and 7 <-> 8. 

After the preparation of the entangled state a waiting 
time of 100 ps was added to let all off-diagonal elements 
decay. After this a pulse of variable length at the partic- 
ular transition was added resulting in the observed Rabi 
oscillations where the initial amplitude is proportional to 
the difference of the corresponding diagonal matrix ele- 
ments. The results were calibrated with respect to the 
previous data obtained for the Boltzmann and pseudo 
pure state. By applying the appropriate numerical anal- 
ysis we extracted for the pp^- state the following values 
for the effective rotation angles aj~ = 86.6°, a} 2 = 86.8° 
and a± s = 88.7°. In this analysis we included the mean 
deviation (e^l = 0.23 which is legitimate since ep enters 
the equations only in even powers. Since all three mea- 
surements are equally sensitive to experimental errors we 
have taken their mean value 

ct x = \ (a- + a{ 2 + a\ s ) = 87.4° . (43) 

By combining the previous results for the initial pseudo 
pure state pro from Eq.(??) with the values of a\ and \ep\ 
we obtain the complete set of the diagonal elements r* 
of the density matrix pp 7 ^- as 

rf = 1-0.01, 0.47, 0, 0.19, i, 0.31, 0.52, 0.02| . 

(44) 

The diagonal elements of the corresponding sublevel 
density matrix p^,- are printed in bold type. Their ex- 
perimental error are estimated to be equal or smaller than 
±0.05. The values are therefore in good agreement with 
the ideal density matrix in Eq. (|32|l . Elements orig- 
inating from the initial Boltzmann density matrix are 
unaffected by the pulse sequence are presented by bare 
numbers or fractions. A similar analysis was made for 
the diagonal elements of the other Bell states. We re- 
mark that this analysis is affected by the decay time of 
the diagonal terms during the waiting time of 100 ps in- 
troduced between preparation and detection. In order 
to estimate its influence we measured this decay time to 
2.60 ms which results in a decay by only 4%. We have 
not corrected for this decay but considered it to be within 
the experimental error. 

3. The off-diagonal part 

Because we have experimentally determined already 
both pulse angles of the sequence for preparing the en- 
tangled states, namely ot\ and f3\ = 7r + e^, the uni- 



tary transformations for the preparation of entanglement 
are known and could in principle be applied to calculate 
the corresponding off-diagonal elements of the density 
matrix. In order to experimentally determine the off- 
diagonal elements we have, however, applied an alterna- 
tive tomography sequence to determine the actual values 
of the off-diagonal elements with respect to the diagonal 
elements. 

In the tomography part of the sequence in Fig. [3] we 
have incremented the rotation angle a (Rabi oscillation) 
at the 7 <-> 8 transition for different settings of the phase 
<j)i and a constant value </>2- After Fourier transform of 
the Rabi oscillations spectra result where their real parts 
are displayed in Fig. El for different values of <fii. Note, 
however, that the signal now depends on diagonal as well 
as on off-diagonal elements of the density matrix. 




200 400 600-0-5 

frequency (kHz) 



FIG. 6: A series of Fourier transformed Rabi oscillations at 
the 7^8 transition of the p^ 9 ~ state for different values of 
<t>i. 

Except for an additional (f> 2 = 90° phase rotation, in 
order to be sensitive to "sin a" -terms, the amplitudes 
were evaluated in the same way as discussed for the Rabi 
oscillations. 

A simple calculation leads to the following expression 

S^- = B + A\ cos a. — A 2 sin + 
A 3 sinasin(20i — <j> 2 ) — A4 sinacos(3</>i — 4> 2 ) ■ (45) 

for the detector signal of this type of tomography. Due 
to the appropriate phase setting combined with select- 
ing the relevant terms of corresponding spectrum after a 
complex Fourier transform allows us extract only terms 
with "sin a" i.e. those two terms with amplitudes ^3 and 
A4 in Eq. 14511 in the data shown in Fig. [7| By using 
the ratio A3/A4 from the phase interferogram in Fig. [5] 
(c) we were able to determine the amplitude A4 from a 
fit to the data shown in Fig. [7| 

The most important matrix element of the entangled 
state pp 7 ^- , namely the matrix element r 2 -j — r-j 2 enters 
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FIG. 7: Amplitudes variation of the Rabi spectra of Fig. [5] 
versus cj>i. 



,_l I L 

1.0. 




-1.0. 



^ 1 1 1 1 

0.25 0.50 0.75 1.00 

time (us) 

FIG. 8: Decay of the off-diagonal element r^i for state 
Pp 7 e _ modulated by double resonance TPPI (time propor- 
tional phase incrementation) with frequencies v\ — 8.0 MHz, 
V2 = 4.0 MHz. The decay is modulated by the characteristic 
frequency 3^i — v% . 



the A± parameter in the following way 
A A = -^r 27 (cos (^)) 3 with ep = ep = ±0.23. (46) 

From the knowledge off all other parameters we obtain 
the value. 

r 2 7 = r 72 = -0.31 ± 0.04 . (47) 

This is only 0.64% of the theoretically expected value 
since from the experimentally known values of the angles 
Pi and a.\ of the preparation sequence and the pseudo 
pure density matrix one would expect = r 7 \ = —0.49. 
The reason for the reduced value is the decoherence of 
the entangled state during and after preparation and the 
delay (including finite pulse width) up to the tomography. 

4- Decoherence 



By analyzing the decay function in Fig. [SJwe were able 
to determine the time constant of the decoherence of the 
entangled state to 

= (208 ± 10) ns. (48) 

Since the process of decoherence is not the subject of 
this contribution we only remark that it is dominated by 
an inhomogeneous distribution of the ESR transitions. 
In order to reconstruct the values of r 27 — rj 2 right after 
their preparation we need to evaluate the effective delay 
time of the tomography sequence. It is comprised of finite 
pulse width as well as actual delay r due to technical con- 
straints. Since the ESR pulses had to be soft in order to 
avoid cross talk with the other ESR line the correspond- 
ing pulse width was 88 ns which is not short with respect 
to . By considering both, the finite pulse width and 
the delay time r we could reconstruct the initial values 
of r 2 7 = r 72 to 



Decoherence reduces the off-diagonal values of the den- 
sity matrix since the pulse sequence for the preparation 
and tomography requires finite pulse lengths and a mini- 
mum delay r between the excitation and the tomography 
sequence. 

We have measured this decay for the VP - state by in- 
crementing the delay time between the preparation se- 
quence and the tomography sequence. The decoherence 
decay shown in Fig. [S] is modulated by the character- 
istic frequency 3^i — v 2 for the pp^_ state imposed by 
the time proportional phase incrementation (TPPI) pro- 
cedure [&2> Lij] • Our type of double resonance TPPI was 
implemented by combined phase increments of <f>\ and 
<t>2 in the tomography sequence. In order to purify the 
decay function from artifacts corresponding to the weak 
lines in Fig. [3] (c) we have added measurements with y 
and — y preparation pulses (phase cycling). As a con- 
sequence the decay function is solely modulated by the 
frequency 3^i — v% characteristic for the entangled state 
\f r_ . Other coherent excitations would appear with other 
phase incrementation frequencies. 



' 27 



initial 



' 72 



= -0.42 



litial 



(49) 



This gives a reasonably high degree of entanglement 
right after creation. 

5. The complete density matrix 

The summary of the density matrix tomography pre- 
sented here results in the complete density matrix for all 
four Bell states. Fig. [5] shows a graphical representation 
of the complete density matrix of pp^ + together with 
complete initial density matrix ppio- 

The elements of the two qubit submatrix are enhanced 
by a darker representation, whereas the other matrix ele- 
ments are presented lighter. Note that the initial matrix 
ppio already contains diagonal elements in the middle 
area which are not relevant for the subsystem we are in- 
terested in. They are only slightly affected by the entan- 
glement sequence. The errors of the diagonal elements 
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FIG. 9: Top: Pseudo pure initial density matrix ppio- Bot- 
tom: Density matrix pp 7 s+ 



of jOp^, + were estimated to be smaller or equal to ±0.05. 
The errors of the nominally zero off-diagonal elements 
are typically smaller than ±0.06. The nominally zero 
off-diagonal elements in the two qubit submatrix have 
errors on the order of ±0.01. Some of the density matrix 
elements are unaffected by the preparation sequence are 
identically zero. 

Similar results were obtained for p^ S/+ and with pp\ 
as initial density matrix for the p]^± Bell states. This 
leads us finally to the following two qubit submatrices 
corresponding to the two qubit Bell states. 



"72 



Off-diagonal elements at detection time: r^. = 
0.31 ± 0.04. Quoted values were reconstructed at prepa- 
ration time due to the measured decoherence time T 9 + 
198 ± 10 ns. 




-0.44 

0.00 0.00 

0.00 ±0.52 



(52) 



Off-diagonal elements at detection time: rf 8 = 
— 0.33±0.04. Quoted values were reconstructed at prepa- 
ration time due to the measured decoherence time T$- = 
210 ± 10 ns. 



r 81 



-0.48 


0.00 
V +0.45 





-0.00 





0.00 ±0.45 



+0.01 0.00 

0.00 ±0.52 



(53) 



Off-diagonal elements at detection time: rf 8 + = 
0.34 ± 0.04. Quoted values were reconstructed at prepa- 
ration time due to the measured decoherence time T$- = 
213 ± 10 ns. 

In order to quantify the accuracy of the experimen- 
tal data with respect to the theoretical expectations we 
define the fidelity as a mean square deviation as 



31 



F P = 1- 



Tr{(p exp - pth) 2 } 
Tr{ P 2 th } 



(54) 



which equals 1 in the ideal case and obeys the relation 
< F p < 1 . The fidelities of all four 8x8 density matrices 
corresponding to the four Bell states are summarized in 
Table I. 











$+ 


$" 


and with ppn 




0.97 


0.97 


0.98 


0.98 



TABLE I: Calculated fidelities of the experimentally deter- 
mined 8x8 density matrices Fp p according to Eq. 1541 . 



-0.01 

±0.47 -0.42 0.00 

-0.42 ±0.52 0.00 

0.00 0.00 0.02 



(50) 



Off-diagonal elements at detection time: 



"27 



' 72 



— 0.31±0.04. Quoted values were reconstructed at prepa- 
ration time due to the measured decoherence time = 
208 ± 10 ns. 



-0.01 

±0.46 ±0.42 0.00 

±0.42 ±0.53 0.00 

0.00 0.00 ±0.02 



(51) 



V. OUTLOOK 

The preparation of the entangled states presented here 
is based on the pseudo pure concept of ensemble NMR 
quantum computing. It has been shown theoretically 
that the corresponding mixed state density matrices are 
separable and do not represent quantum entanglement in 
the strictest sense (43, |43| • Even if we would deal with 
a single 15 N@Cqo molecule the pseudo entangled states 
presented here would be mixed states at the temperature 
of T = 50 K and microwave frequency v — 9.5 GHz ap- 
plied here. However, in case we would start out from a 
reasonably pure state and perform the same preparation 
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and tomography scenario as presented here we would in- 
deed reach a quantum entangled state. No other unitary 
transformations as presented here would be required and 
in fact the same signatures of the tomography like the 
characteristic phase dependence would be observed. In 
this sense we consider the experiments presented here as 
precursors of the corresponding quantum experiments. 
In the following we try to estimate which experimental 
parameters are required to reach such a quantum state. 
We start from the Boltzmann density matrix according 
to Ea. l|2U|l where the inverse temperature j3 = h/(k-g,T) 
and the Larmor frequency of the electron spin ll>s = 2irvs 
play the dominant role. We neglect any Boltzmann polar- 
ization of the nuclear spins because of their low Larmor 
frequency. In order to produce a spin alignment as initial 
state we first apply a (3q = tt pulse at the electron spin 
transitions with mj = —1/2. On this initial state we 
apply the unitary sequence according to Eq. (J2SJ. This 
leads to a density matrix which supposedly represents a 
quantum entangled state for certain values of tempera- 
ture T and frequency wj. All eight eigenvalues of ps are 
positive at any frequency lus and temperature T . 

In order to apply the positive partial transpose (PPT) 
criterion of Peres and Horodecki 0, [5{| if a mixed state 
density matrix is separable or not we need to perform the 
partial transpose on one of the spins. If all eigenvalues are 
still positive after partial transpose, the density matrix 
is separable. If at least one of the eigenvalues becomes 
negative there is some degree of entanglement. After 
performing the partial transpose on the nuclear spin / 
we obtain a zero crossing of one of the eigenvalues under 
the condition 

i 

e ^ s >(s + 2V2) 3 or pu s > - log (3 + 2V2) . (55) 
The value for the quantum critical temperature below 



which the quantum regime is reached can be expressed 
as 



1 Q = 1 7 FT- °°J 

q fc B ilog(3 + 2\/2) 

For a high field ESR spectrometer operating at 95 GHz 
one obtains T q — 7.73 K. This is well in reach of cur- 
rent experimental setups. It demonstrates that with- 
out the usual dynamical nuclear polarization (DNP)of 
nuclear spins quantum states can be created and corre- 
spondingly quantum algorithms can be performed. Ex- 
periments along these lines are planned. 



VI. SUMMARY 

We have shown in detail experimental results on how 
to prepare pseudo entangled states in a confined spin 
system, namely 15 N@Cqo consisting of an electron spin 
S = 3/2 and nuclear spin 7 = 1/2. All four Bell states of 
a two qubit system were produced. Moreover, we have 
performed several variants of spin density matrix tomog- 
raphy in order to verify the degree of entanglement. By 
estimating the quantum critical temperature T q — 7.73 K 
for an ESR frequency of 95 GHz we could show that true 
quantum states are in reach with current technology. 
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